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Abstract 

The flux of hypothetical "hidden photons" from the Sun is computed under the assump- 
tion that they interact with normal matter only through kinetic mixing with the ordinary 
standard model photon. Requiring that the exotic luminosity is smaller than the standard 
photon luminosity provides limits for the mixing parameter down to x ^ 10^^*^, depending 
on the hidden photon mass. Furthermore, it is pointed out that helioscopes looking for solar 
axions are also very sensitive to hidden photons. The recent results of the CAST collabora- 
tion are used to further constrain the mixing parameter x low masses (m-y/ < 1 eV) where 
the luminosity bound is weaker. In this regime the solar hidden photon flux has a sizable 
contribution of longitudinally polarized hidden photons of low energy which are invisible for 
current helioscopes. 



1 Introduction 



The standard model (SM) provides an accurate description of particle physics below the elec- 
troweak scale but it is generally though not to be valid up to arbitrary energies. Extensions 
of this celebrated scheme, invoked to cure diseases like the strong CP, hierarchy or flavor 
problems, often involve higher gauge symmetries and further matter content. Moreover, 
string theory is a preferred candidate for the unification of quantum mechanics and general 
relativity where additional gauge and matter fields are assured. At low energies, some of 
these new fields can arrange into a "hidden sector" if only very massive particles (or gravity) 
mediate interactions between them and the SM "visible sector" . 

Of course, depending on the scalar content of the theory, gauge symmetries can either 
be spontaneously broken or remain exact. Then, the corresponding "hidden" bosons could 
have in principle an arbitrary mass. If this is small enough, these hidden bosons can have a 
very rich phenomenology at present affordable energy scales. 

The simplest case concerns just a novel f/(l)h symmetry and its corresponding gauge 
boson, henceforth called "hidden" photon. The interplay between this hidden photon and 
the SM photon modifies the predictions of quantum electrodynamics [1], often claimed to 
be the most accurate of all physical theories so far, thus constraining the hidden photon 
parameters. We can turn this argument in the opposite direction: the constraints on hidden 
photon parameters give us information about how accurate is the QED description of nature 
at low energies. 

A number of laboratory experiments has been devoted to the search of hidden photons, 
the resulting bounds being strongly dependent on the hidden photon mass. For masses 
corresponding to macroscopic length scales, experiments testing the Coulomb law |2l |3] set 
strong constraints on hidden photons, but still they could be largely improved by experiments 
dealing with high quality microwave cavities [1]. In the microscopic range, laser experiments 
are also becoming very powerful probes of hidden sector particles [3 El [71 El El [iDl ttH US] . 
At atomic distances, comparison of the Rydberg constant for different atomic levels gives 
interesting but weak bounds [T3l [Ti] . Finally, particle colliders extend the mass range until 
typical electroweak scales [ISl [IHl ttH [IHl [IH] • 

On top of that, the evolution of stars turns out to be the most sensitive "laboratory" 
to study properties of novel low mass weakly interacting particles [201 EI] • Even with tiny 
couplings to electrons and protons, they might be still copiously created in the interior 
of hot and dense stars. Because of their weak interactions they might abandon the star 
without further scattering, accelerating the consumption of nuclear fuel, and therefore the 
stellar evolution [22l [23]. Our present observational data on stellar evolution can strongly 
constrain this novel luminosity, although the bounds can be relaxed in some concrete models 
[21 [251 [261 [271128]. 

Interestingly enough, in this case the Sun itself could be a copious emitter of weakly 
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interacting particles, that could eventually be detected at Earth inside a sensitive detector 
[25] (as this is actually the case with neutrinos). Several of these so-called "helioscopes" 
[30] have been built with the aim of detecting solar axions [311 ESI |33] and remarkably, the 
CAST collaboration has even recently surpassed the sensitivity of the energy loss arguments 
for the axion coupling to two photons [3l]. As we will see, helioscopes can also detect hidden 
photons so the CAST limits can also be used to constrain the solar hidden photon flux. 

The energy loss argument and helioscope bounds for hidden photons were already studied 
in [55| EHj- However, in this paper the author does not consider either the possibility of a 
resonant production, which can enhance enormously the hidden photon flux, nor the emission 
of transversely polarized hidden photons. 

In this paper we compute the energy loss bounds using the latest solar data [37] and 
derive the CAST helioscope bounds. In particular, it is shown that the new CAST results 
are extremely sensitive to solar hidden photons, providing the strongest constraint of their 
existence in the mass range m^/ ~ 0.01 — 1 eV. Moreover, accounting for the resonant 
production improves the energy loss bounds in [36] roughly up to 1 order of magnitude. 

The paper is organized as follows: In Secj2] the hidden photon solar emission is derived 
while the principles of the CAST helioscope detection are reviewed in Sec|3} We compute 
the energy loss and CAST bounds in Sec|4]and finally the conclusions are presented. 



2 Hidden photon production in the Sun 

At low energies, the dynamics of the two photon system can be described by means of the 
following Lagrangian, 

which is the sum of the Maxwell Lagrangian for the standard model photon field with its 
corresponding source (the electric current of electrons and protons), the Proca Lagrangian 
for the massive hidden photon 5^ and a gauge invariant mixing term. This term has been 
discussed in some detail in the literature in a variety of contexts [3H1 ESI SOI ISl 1121 SSI 
HH HHl Hni ST]. Although X is a priori a free parameter that could be even of order 1, the 
experimental evidence indicates that it should be much smaller. A natural explanation for 
this is that while x could be zero because of symmetry reasons at a high energy scale, the 
integration of high energy quantum fluctuations will end up inevitable in a small nonzero 
value in the low energy theory. Typical values comprised in the literature range from 10^^^ 
to 10~^ [H]. The lack of a mass mixing term is of course due to the protection of the f/(l)em 
symmetry, which in principle we would like to preserve. We do not have to write a tree level 
coupling B^j^ because it can be eliminated (if small) by a harmless redefinition of the photon 
field and therefore absorbed in x- This could not be done if the B^ field couples with different 
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strength to electrons and protons. Such a situation, however, is severely constrained from 
experiments testing the neutrality of matter [IH] and arguments concerning the existence of 
leptonic and baryonic forces [IHj and therefore we will not consider it. 

The presence of a mixing term simply states the fact that the and B^^ fields are 
non-orthogonal. This can be easily seen if we notice that the Feynman rule for the A^^B^^ 
vertex is simply the inverse (non diagonal) massless propagator. In order to gain further 
understanding of the physics to be discussed it is highly recommendable to express the 
Lagrangian in terms of the state that is orthogonal to the photon, and therefore sterile with 
respect to local electromagnetic interactions, defined as 

= 8^" + xA" . (2) 

This invites us to renormalize the electric charge a/1 — x"^^^ ~^ and generates non- 
diagonal mass terms, 

= \m\S^S, - xmt^'S^A^ + ^x'K'^^'^m > (3) 

which indicate that photons oscillate into sterile states during free propagation, like in the 
case of neutral kaons or neutrinos. 

The equations of motion (EOM) in the Lorentz gauge for such a system are, in Fourier 
space, 

{K^g^^ -W^{K)-x^nil^>)A,{K) + xml^^S^{K) = , 

{K^-m%)S^{K) + xml^,A^{K) = 0, (4) 

where K is the four momentum (from now on I will skip writing it) and I have included 
the effects of the photon interactions with the medium via the polarization tensor Tl^^ , 
which has a clear physical interpretation as the photon self energy in the medium]^ In a 
homogeneous, isotropic and CP conserving medium it has only two independent components 
,7rT and ttl, corresponding to two transverse (k ■ A = 0) and one longitudinal (Langmuir 
waves with k x A = 0) electromagnetic excitations, which in a plasma are called plasmons. 
Taking advantage of the isotropy, we chose to focus on propagation along the z-direction 
K = {uj; 0, 0, k) and we can write a simple set of orthogonal polarization vector^ as 

e^^(0;l,0,0) ; e^^(0;0,l,0) ; e^^=L=(A;; 0, 0, ^) , (5) 

satisfying e^e^^ = —Sab and define 

U>^'^(K) = -J2<C^aiK) . (6) 



^In this discussion I follow the exposition and notation of [20] , 

^We will not consider space-like plasmons in this paper so the denominators will cause no trouble. 
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Thus we can write the EOM for the different components Aa = ^a^f^ (idem for Sa) as 



The evolution of T and L plasmons is decoupled. 

The lowest order contribution to 11'^^ in a plasma comes from coherent forward scattering 
off the thermal bath of electrons [50l|5T], providing refractive (real) parts for 7rT,L- For typical 
solar plasmon energies, < 10keV<C rrie we can neglect the electron velocity dispersion in our 
calculations getting [51] 



Re{7rT} = m| ~ cul ; Re{7rL} = ttil ~ - A;^ , (8) 
with the plasma frequency up given by 

Up ~ Tie , (9) 

me 

where He is the electron number density. 

Leaving aside the effects of the hidden sector photons, the solutions from the equations of 
motion eq. ([T]) tell us that T-plasmons behave as massive particles with an "effective mass" " 
given by up but L-plasmons oscillate at a frequency u up, almost independent of the 
wave number k (the group velocity is suppressed by a small factor ~ 3T/me). Then, while 
T-plasmons are always time-like, L-plasmons can also be light or space-like. Explicitly one 
finds 

ST 

lj^c:^ujI, + e (T) ; ^ 4 + —k^ (L) , (10) 

TTle 

Where I have kept the lowest order terms in k. 

In order to get the lowest nonzero contributions to the imaginary parts of Ht^l we need 
to include photon absorption and dispersion. Following Weldon [52], the imaginary part of 
the photon self energy is proportional to the difference of the photon absorption and photon 
production probabilities (F^'^) by means of 

Im TTa ^ -iuTa = -Lu {Tt - FD = -uj{l - e-^)F^ , (11) 

where in the last equality thermodynamic equilibriun:(^ is assumed, implying 

FP = fV? . (12) 



The result eq. (11) reminds on an absorption rate corrected for stimulated emission. 

The dominant source of opacity in the Sun is inverse bremsstrahlung (also called free- 
free absorption) which is very efficient at low and intermediate energies. Compton scattering 



^Even in the absence of thermodynamic equihbrium the result holds because of unitarity |53L I52j . 
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provides a smaller but energy independent contribution (given that uj <^ irif.) which is crucial 
for the highest energies. For transversely polarized photons of energy sufficiently above the 
plasma frequency we get 



where is the number density of ions of charge eZi -which as a good approximation we will 
take to be only Hydrogen and Helium- and ~gffi is the Boltzman averaged Gaunt factoin 
which accounts for the deviations from the classical expression derived by Kramers [57]. 
more accurate calculation should include the effects of free-bound and bound-bound tran- 
sitions, which might acquire some relevance at the relative low temperatures of the solar 



external layers. Note that in eq. (13) we are also assuming that all the particle species 
are completely ionized. Again this approximation will fail at the solar external layers and 
corrections should be included by solving the Saha equation. We will refer to the validity of 
these simplifications later on. 

For ttl the situation is more complicated since formulas for the interactions of longitudinal 
plasmons are difficult to derive or to find in the literature. However it turns out that in the 
most interesting case, where L-plasmons are time-like, the bulk emission is not sensitive 
to the details of Fl. For completeness, however, we can write the simplest contribution, 
Thomson dispersion (derived in Appendix |A]), giving 

„ Svra^ k , 

Ue . (14) 



9meT uo 



The solar model BP05(OP) [37] provides the latest available data on the solar interior 
compatible with helioseismology and neutrino fiuxes. Magnitudes like temperature (T), mass 
density (p) and the mass fraction of the most important atomic elements (Xj) are tabulated 
as a function of the solar radial coordinate r. Assuming that the plasma is locally neutral we 
can compute the electron density as the sum over species of the number of protons through 
rif. = p/niu'^^ZiXi/Ai, with Zi,Ai the atomic and mass number of species i and the 
atomic mass unit. The relevant parameters for this work are plotted in Fig. [T} In particular 
notice that the electron density strongly depends on the position in the solar interior. This 
introduces an implicit position-dependence in the EOM ([T]). 

Now that we have learned about the standard plasmon propagation in the Sun it is time 
to address the effects of mixing with the hidden photon B. Assuming slow spatial variation 
of the electron density, the EOM in eq. ([T]) are diagonalized by the following shift 

2 2 
Aa > Aa + X ] Sa > iSq + X 9 (1^) 

''For a review on the subject see [51] . In the numerical estimates of this paper I have recalculated the 
Gaunt factor using the exact Sommerfeld-Maue formula |55j and a Laguerre-Gaus quadrature as proposed 
in [SB]. 
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Figure 1: Values of the solar parameters relevant for this work plotted as a function of the 
normalized solar radial coordinate r (left) and the plasma frequency up (right) in decimal 
logarithmic scale. From small to large dashing these are the electron density rig (red), 
temperature T (blue), plasma frequency (black), dup/dr (brown). Hydrogen mass fraction 
X (green) and radial coordinate r (pink). Except for X, they are normalized to their largest 
values 6.07 x lO^^ cm'^, 1350 eV, 295.5 eV, 5.8 x 10"^ eV^ m'^ and Rq = 6.96 x 10^ m, 
respectively (taken from [37]). 



as long as the non diagonal elements are much smaller than 1. I will refer to this condition 
as weak mixing (WM) and justify later that it is satisfied in all the relevant cases. 

The states A, S have decoupled evolution and satisfy dispersion relations similar to A 
and B in the absence of kinetic mixing, 

u;'-kl = myliu;,ksJ=m'^, + 0{x') . (16) 

Therefore an originally pure-plasmon state A with energy u produced in a small region 
labeled z = will evolve as a linear combination of the propagating states, 

Aa{t, z) ~ e*-*-^^-4.M„ + , f"",^' ^e'^'-'''''^'~Sa , (17) 

leading in principle io A — S oscillations. However, in the Sun the imaginary part of the 
plasmon dispersion relation is big enough to damp completely the A component after a short 
distance ~ F"^. 

The S component is damped as well since it receives an 0{x^) imaginary contribution 
oc Fa. Once we know what values of x "we can bound we can go back to this point to show 
that this absorption is negligible. 
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Finally we have to take into account that this solution holds locally in the Sun, or at 
least in a region where rif, (and therefore tTo) does not change very much. At every position 



of the Sun we should define Sa{r) inverting eq. (15) and using eqs.(|8])-(14) with the function 
ne(r) given by the solar model. Interestingly enough, as long as Ta is sufficiently large, 
the corresponding Aa{r{z)) component will be quickly absorbed and in practice Aa{t, z) will 
"follow" S'a(r(z)) as it travels out of the Sun. 

Eventually, a plasmon will exit the Sun as a Sa^Ro) state, 

~Sa{Re) = ~Sa{ne{Re) ^ 0) = 5„ , (18) 

with a probability 

4 

Pa^^b. = \{SMa{Re))\' ^ X'^, 2 J.' , . ^. . (19) 



where m^Q = m^(ro) and ci;ra,o = ^a{ro) are the real and imaginary parts of tTq evaluated at 
the production point, at a distance tq from the solar center. 

Let me remark that the hidden photons exiting the Sun will have the same energy than 
their original plasmons, i.e. of order eV up to ~10 keV. As Ba are propagating states in 
vacuum they will move away from the solar surface without "flavor" oscillations. 

The rate of hidden photons exiting the Sun will be then proportional to the plasmon 
generation rate, i.e. the rate at which plasmons are produced from electron and proton 



interactions in the plasma. Such a generation rate can be derived from eqs. (13) and (14) by 
using eq. (11) and eq. (12). We find that this equal^Fa defined in eq. (11) weighted by the 



Bose-Einstein distribution ube = nBEi^,T) = [er — 1] . This is particularly transparent in 
the case of scattering since every photon emitted comes from a thermal photon "absorbed" 
and the number of these photons is given by the Bose-Einstein function. 

We can immediately compute the fiux of hidden photons that arrive at the Earth by 
integrating this rate over the solar model [37] , 

where one has to keep in mind that and Ta depend implicitly on r through ne{r) and 
eqs.(|8|-(14) and on u; = uj{k) given by eq. (10) . i?© is the average Sun-Earth distance 
~ 150 X 10^ m. We discuss separately the emission of T and L hidden photons because the 
different dispersion relations lead to a completely different phenomenology. 

^This is of course related to KirchhofF's law of thermal radiation. 
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2.1 production 



Using the dispersion relation in eq. (10)(T) and multiplying by a factor of 2 to account for 
the 2 different T-polarizations, we get 



d<l> 



T 



dtu 



4vri?| JO 



R 



I 

Aur'^dr — 



2 4 



CT — 1 



(21) 



Given that in the solar model 1 eV < cup < 295 eV, this expression has three clearly 
differentiated regimes that we discuss separately. 



2.1.1 Suppressed production (my <^ 1 eV) 

For m-y' <^ 1 eV we can safely neglect the terms involving my and Ft (typically Ft < cup) 



in the denominator of eq. (|19|) and in the square root. Then, all dependence on the hidden 

A' 



photon parameters, x^^y? factors out of the integral in eq. (21) and we have 



d<^> 2 4 1 /"^® 2,1 Tt o fmy\^ FUu) , , 



dio ^ 47ri?^ Jq tt^ eT — 1 tUp v eV / cm^ s eV 

where the dimensionless function -Fi(u;) is plotted in Fig. |2| The spectrum peaks at low 
energies because of two reasons: on the one hand the bremsstrahlung production decreases 
as a;~^ and on the other hand the production at the solar core, from where the most energetic 
photons are expected, is suppressed by the largest values of uj-p. 

As we will see, this suppression regime is the most interesting for helioscopes. A simple 
analytical formula for Fx can be useful therefore for more delicate future analysis. A fit 
similar to the axion flux in has provided, 

Fi(cj) = 2.7 X 102*^E^2.98g-^ ^ (23) 

where E = u/keV. The fit reproduces the numerical results between 0.5 and 5 keV with 10% 
accuracy. 

The solar luminosity in hidden sector photons in this case is 

= xX' r ^^r'dr\^^^ = 8 x 10^^ C^Y Watt . (24) 

It is interesting to point out that for low energies ~ eV most of the production comes 
from the outer layers of the Sun. In this region, neglecting bound-free and bound-bound 
transitions or the ionization fraction is not completely justified. However, on the one hand 
they contribute relatively less than higher energies to the energy loss, and on the other hand 



the CAST bounds (based on keV energies for which eq. (13) is accurate) are more restrictive. 



Therefore, unless this uncertainties imply a huge increase on the energy loss, which is not 



likely, or an helioscope focuses on these energies, there is no need of correcting eq. (13). 
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WkeV <^eV 

Figure 2: The Fi and F2 functions give the flux of solar transverse and longitudinal S's at 
the Earth for m^i -C 1 eV. Notice the different energy scales, only eV L-hidden photons are 
emitted while the spectrum of T-modes extends to X-ray energies, although considerably 
suppressed. See the text for details. 



2.1.2 Resonant production (leV< my < 295 eV) 

In this mass interval, there is always a small region in the Sun where the hidden photon 
emission is so intensely amplified that outshines the emission from the rest. This is the 
region where the plasma frequency is tuned to the hidden photon mass uj-p = m^/, and 
correspondingly the emission probability is 

PaZb = . (25) 

The factor cuFt/cUp, plotted in Fig. [sj enhances the mixing and therefore the probability, 
but if it is too small can invalidate the WM condition. We will find that the energy loss 
argument imposes values for x smaller than 10~^ in this mass range (See Fig. [6]) so in any 
case the WM condition is always satisfied. 



If we neglect the small r-dependence of T and Ft in eq. (21) the half- width of this 



resonance can be easily calculated from eq. ( 19 ) giving 

Acu| = 2wFt . (26) 

Therefore the resonant emission will take place in a tiny shell of the solar interior, of size 
Ar ^ /S.ujp{duj'^ / dr)^^ . The resultant flux at the Earth can be approximated by 
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Figure 3: The value of ujT/uj'^ controls the enhancement of the probability of emission of 
hidden photons in resonant conditions. The shadowed region contains the values for the 
transverse modes in the whole solar model, with the boundary curves for the Solar center 
(up) and surface (down). The short line is a lower bound {m^i = up) for the longitudinal 
modes for which the relation between the energy and the position at the Sun is fixed. See 
the text for details. 



As a remarkable fact, the dependence on Ft cancels out. Any improvement in the derivation 
of Ft will not change this result. However, even coming from a tiny shell of the Sun with 
almost constant temperature T = T(r (tup)), the spectrum of hidden photons does not have 
a perfect thermal shape; the power of the energy is roughly u, not u"^. 



Finally the energy loss in hidden photons can be easily calculated. Using u!j> 
I find 



duj 



-ujduo 



16C(3) r^T^ 



dr 



4 2 



{2i 



Let me remark that while for m^i ^ 1 eV the m^/-dependence of the hidden photon emis- 



sion factors out, here in eqs.(|27|) and (|28|) it is implicitly assumed that all the quantities 

These 



my. 



T, r, (x'P, Ft, (icjp/fir are evaluated at the point of the solar model where tup 
dependences can be read in Fig. [T]and will help us to understand the bounds of Sec. |4| 



2.1.3 Unsuppressed production (m^/ > 295 eV) 



In this case the emission probability of eq. (19) is simply x^. The flux in eq. (21) is therefore 
independent of m^i if the energy uj is high enough such that the threshold corrections like 
the square root are small. In this case we find 



rf<l>T 
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UJ 



cm^ s eV 



(29) 
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Figure 4: The function G gives the flux of m^i ^ 295 eV hidden photons from the Sun. See 
the text for details. 



where G{(jj) is a dimensionless function plotted in Fig. |4j In general in this regime most 
of the production comes from the solar center, where Ft and T are bigger and therefore 
also the possible energies uj. For energies above ~ 5 keV Compton scattering dominates 
over bremsstrahlung and thus we expect an almost Planckian spectrum with a temperature 
~ T0 = 1.35 keV, with exponential suppression for m^i 3> Tq. 



2.2 5l production 



Longitudinal plasmons can also resonantly convert into hidden photons. Moreover, since in 
this case ~ tUp — /c^ the resonant conversion will not only happen at a solar shell where 



but at every place in the Sun where uj-p > m^i since mi^ = mt is always satisfled 



by a certain value of k. 



As the energy ~ cup is almost unrelated with k by eq. (10)(L) but r = r{ujp), it is 



convenient to approximate the momentum integral in eq. (20) by the value at the resonance 



times the half-width Ak (in analogy with eq. (26) Afc^ will be 2ci;Fl) and use 



dr 



dujp 
dr 



-1 



2ujduj 



dujp 
dr 



-1 



(30) 



to get 



doj 



Ra 



TT 



2 d.^'i' 1^ 
dr T 



(31) 



which looks exactly the same formula than for the resonant flux of T-hidden photons eq. (27) 
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but here r and dup/dr are evaluated at the point where u-p = uo. This hmits the Bl flux to 
energies leV ^ < 295 eV, the range of the plasma frequency in the solar model. 



For small m^/, again m^x^ factors out of the flux and, in analogy with eq. (22), we can 
define a dimensionless function F2{uj) as ((i$L/c?^^)(cm^eVs)(my/eV)~^X~^ which is plotted 
m Fig. [2j In this mass regime the 5l luminosity is 

WT{my < 1 eV) = 1.4 x 10^^ (^j Watt , (32) 



which is only slightly smaller than than the luminosity in eq. (24). 



Then, while resonant T-hidden photon production proceeds in a tiny shell that emits at 
all energies, L-hidden photons are resonantly produced in a sphere of radius r = r^u-p = my) 
and the position inside this sphere determines the unique resonant energy ut = ujp{r). 

Unfortunately, beyond m^/ = lupq ~ 295 eV the resonant production is not possible and 



the emission of -Bl's drops drastically. Note from eq. (10) that the L-plasmon energy cannot 
exceed cup by much without entering the region where Landau damping is strong, namely 
k ~ \/'nT'e/Tujp. For such large values of k all our approach has to be revised. However, at 
the view of the bounds that can be extracted with the more powerful resonant emission this 
seems not to be a fruitful business. 

It is easy to show that the WM condition in this case is also satisfied. For my > 1 eV we 
will require (at least) x < 10"® and Fig. [s] shows than at most an enhancement of lO'* can be 
expected. For smaller values of m^/ these quantities both scale with my so the conclusion 
remains unchanged. 



3 Helioscope detection 

Having addressed the calculation of the hidden photon fiuxes at the Earth, in this section 
I point out that existing axion helioscopes [29l |30] like CAST |33l [31] at CERN are indeed 
capable of detecting these exotic photons. 

The set up of a typical axion helioscope is depicted in Fig. |5} It consists in a (preferably) 
long cavity pointing towards the Sun where a strong magnetic field is maintained. The cavity 
is strongly sealed and has a powerful, low background X-ray detector at the end. 

If axions are emitted from the Sun [221 EH] , they easily pass through the shielding and a 
few of them can be coherently reconverted into X-rays. This conversion can be understood 
in terms of axion-photon mixing [59], where the magnetic field acts as the mixing agent, 
providing the required angular momentum (axions are spin-0 particles). 

In contrast, hidden photons do not need a magnetic field to satisfy angular momentum 
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Sun 




Figure 5: Schematics of an helioscope experiment looking for hidden photons. Hidden bosons 
are produced in the Sun's interior from X-ray plasmon conversion and get out almost freely. 
Only the sterile component (5*) of B will traverse the helioscope external shielding, leading 
to the possibility of reconversion into a detectable photon hj S — A oscillations. 



conservation, they just mix naturally with photons regardless of the presence of the magnetic 
field, given that x ^i-nd m^/ are both non zero. 



The small mixing shift in eq. ( 15 ) turns out to be again a correct approximation for the 
treatment of the photon-hidden photon system both in the Sun-Earth travel and inside the 
CAST oscillation region. This is because CAST operates either in high vacuum where we 
can neglect TTa (and the WM condition is an exceedingly good approximation), or filled with 
gas to force resonant conversion, and therefore producing also a non-zero Ta that cuts off 



the divergence of the denominators in eq. (15) at reasonable values. The interstellar medium 
can be also treated perfect vacuum. 

Note that essentially all the hidden photons emitted from the Sun are propagation states 
oc Ba so they do not suffer oscillations. However, these states will be projected into their 
sterile component S when traversing the helioscope shielding. Such a projection decreases 
the flux only in a small O(x^) factor and can be neglectecj^ An initial state Sa traveling 
through the conversion region (of length L) will oscillate into a detectable photon with a 
probability given by 

Ps.^A. = \{Aa\Sa{L))\' = \{A^{L)\Sa)\' = (33) 

^'"^^ ^l + e-r''^-2e-^cosAp„L) , (34) 



(m2-m2,)2 + (^r,)2 



where Apa = y^u^ — — yuj"^ — to^, is the difference in wavenumbers of the photon and 
hidden photons. In vacuum = F = and L-plasmons cannot been excited, so we recover 



^Effects of the atmosphere and further barriers are also at the ©(x^) level and therefore are also unim- 
portant. 
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the well-known expression 

P^^^A^ =4x'sin2^ , (35) 



that in the limit uj ^ my leads to 



m.L 



^Sx^Ax=4x'sin2-^ , (36) 
4a; 



and, for my ^ ujL ^, 



An interesting situation arises when the hidden photon mass is so large that the argument 



in the sinus of eq. (36) is much larger than one. If the energy resolution of the detector is 



such that integrates photons with energy between uj and uj + Aa; then it effectively integrates 
several oscillations when Acj/cj ^ AniLO j{vp?^,V). In the limiting case, the sinus is effectively 
averaged to 1/2, if the dependence of the flux with the energy is reasonably small. Note that 
in this case the conversion probability will be effectively independent of the oscillation path. 

If oscillations take place in a medium, even longitudinal excitations could appear. How- 
ever the detectors used for helioscopes like CAST are designed for the detection of trans- 
versely polarized photons. For the T-modes a resonant detection is possible if the oscillation 
volume is filled with a small amount of gas such that m^(= cUp) = (above the energy of 
the highest atomic resonance the dispersion relation of photons in gas is essentially the same 
than for a plasma). Assuming V ^ 1, the conversion probability is again independent of 
the absorption coefficient, 

A formula for the range of validity of the WM approximation can be easily provided. For 
keV photons above the atomic resonances of the gas the Thomson scattering will be the main 
source of absorption. We find then that x^p/ {^'^t) = x3'^e/(2aco') ~ 10^x(keV/ci;) <^ 1 
should hold. 

In our bounds we should use the values for the most recent and sensitive experiment. The 
CAST helioscope at CERN has recently published results |3l] from a search of solar axions 
in a decommissioned LHC magnet of 10 m length in an energy window 0.5 keV< u; < 15 keV. 
They found no signal over a subtracted background of 1 ~ 4 x 10~^ counts cm~^s~^keV~^ 
depending on the energy. In order to be conservative I will use a slightly bigger value 

$CAST < 10~^cm~2 s-^keV~^ , for 0.5 keV < to < 15 keV , (39) 

in the bounds of this paper. This limit arises from the difference of background photons 
when the helioscope points towards the Sun and when it does not (axions need the magnet 
to convert into photons). If the hidden photon mass is big enough such that the conversion 
probability is independent of the length this procedure will subtract the possible signal as 



well! Therefore, for masses larger than ^ySiruj'^ /{AloL) one should use the pure background 
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Figure 6: Upper limits on the kinetic mixing parameter vs. the mass of the hidden photon 
from the CAST hehoscope and the Solar lifetime argument. See the text for details. 



counts to estimate the bounds. However, one of the 3 detectors of the CAST experiment 
(the CCD camera) has a X-ray focusing device that in any case increases the flux of photons 
coming along the magnet direction with respect to those coming from other directions in 
a sizable amount, making the subtraction harmless. As the exclusion bound derived with 
this detector is stronger than with the two others (the TPC and the Micromegas), it seems 



conservative to still use eq. (39) to limit the flux of photons from hidden-sector conversion. 
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4 Bounds and discussion 



The energy loss in terms of novel particle species can have dramatical consequences for 
stellar evolution. Hidden photons seem to behave like invisible axions [60l EH [29] , wherever 
they are produced they leave the Sun without further relevant interactions. Although in 
principle such a non-standard energy loss could be easily accommodated in a "present day" 
solar model by increasing the central temperature a bit over the standard value (the nuclear 
reaction rates depend strongly on the temperature) the nuclear fuel would be consumed 
faster and the star would enter sooner in the next stage of stellar evolution [23]. 

Theoretical [23] and numerical studies [62] incorporating axion loses showed that no 
present day solar model can be constructed if the exotic luminosity is bigger than the actual 
solar luminosity Lq. Under the hypothesis that hidden-photon losses have similar effects 
and taking Lq = 3.83 x 10^^ Watt we can obtain the bounds in Fig. [6] labeled T-Lifetime 
and L-Lifetime. 



The T-Lifetime bound has been obtained by integrating numerically the general expres- 



sion eq. (22) for Bt, except for the region of resonant production (1 eV< m^/ < 295 eV) 



for which eq. (28) has been used. The L-Lifetime bound comes from numerically integrat- 



ing the energy of the flux in eq. (31). The bound should regard the sum of both T and L 



contributions but in practice the L contribution does not change the limit appreciably. 

Below my ~ 1 eV the bound is a straight line of slope —2 because in this case the hidden 
photon flux is just proportional to x^'^y- the range 1 eV < my < 295 eV the resonant 
production dominates and improves the bound down to x ~ 10^^"^ where the resonance takes 
place in the solar core. The change of slope of the T-Lifetime bound around m^/ ~ 10 eV can 
be traced back to the temperature dependence of eq. (24), as can be seen in Fig. [l] since in 
this regime up corresponds to m^/. There we can also see that the inflection point at ~ 290 
eV it is due to the decrease of the volume of emission through the factor of eq. ([24]). 

Beyond 295 eV the T-Lifetime bound sharply worsens until it reaches the bound due 
the solar bulk emission. I 
photons with enough energy. 



to the solar bulk emission. For my ^ 10 keV the bound vanishes because of the lack of 



Let us now turn our attention to the CAST bounds coming from the non-observation of 
X-ray photons in the energy window 0.5 — 15 keV. We have already discussed that CAST 
is not sensitive to -Bl's so we restrict the discussion to the non-observation of transverse 
excitations. 



The hidden photon flux is a decreasing function of the energy since it comes essentially 
from bremsstrahlung (see also eq. (27) and Figs. |2]and|4]) so we can restrict the bounds to 
the most interesting (yet reasonably sized) interval, which is 1 keV above my. 



Integrating the flux in eq. (21), or eq. (27) when resonance dominates, times the 
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vacuum conversion probability eq. ([35j) in the mentioned energy interval and using eq. (39) 
I get the bound labeled "Helioscope (CAST)" in Fig. M 



The onset of the oscillation regime in the CAST detector is clearly seen at m. 



■7 



^4 X lkeV/10 m ~ 0.01 eV. Below this mass, S — A oscillations are never complete and 



eq. (37) holds. The CAST signal is then proportional to x mL, and therefore the slope of 



the Helioscope line is again —2. Soon above this mass, the oscillation length in CAST is so 



small that the squared sine of eq. (35) gives a factor of 1/2 when we average over energy. 



In this regime, but below 1 eV, the CAST signal is proportional to x™y ^ind therefore 



the Helioscope line has slope 1. In order to mantain this slope down to 10"*^ eV the CAST 
oscillation length should be increased in a factor 10^ (or the energies detected lowered in a 
similar factor), however the region of improvement is already excluded by laser experiments 

m. 



Assuming that the same bound eq. (39) holds even with gas filling the oscillation region, 
and varying the plasma frequency of the gas in the reasonable range 0.01 — 1 eV, the bound 
labeled "Gas" could be achieved. While it is not clear if this procedure is realistic for higher 
masses. Fig. |6] shows that it would be extremely interesting. 

In the mass range ~ 1 — 10^ eV, the CAST bound lies in a region ruled out by the energy 
loss argument. Since the emission has been calculated assuming an unperturbed solar model 
we find that the CAST bound is not consistent in this mass range. The observation of a 
hidden photon flux corresponding to this region would imply drastic changes in our current 
picture of solar structure. 



Note that the curve scales with the fourth power of the CAST limit eq. (39) so a huge 
improvement in sensitivity would be needed to beat the energy-loss bound in this range. 
Observe however, that there is no change of slope in the Helioscope line at ;U ~ 1 eV where 
the T-Lifetime starts to be dominated by resonant emission. The reason is that the resonant 
production starts at cup ~ 1 eV, close to the solar surface where the temperature is much 
smaller than the CAST energy window. Therefore, some hope of improvement in the mass 
range 1 — 10 eV relies in lowering the lower threshold until (9(eV) energies. 

The top of the energy window, 15 keV, limits the range of hidden photon masses that can 
be testable with CAST. The dashed vertical line shows this limitation. Still I have continued 



the "Helioscope" line up to somehow higher masses by assuming the same sensibility eq. (39) 
and no threshold. Interestingly enough, this shows that for m^/ > 25 keV the CAST bound 
would be again consistent with, and more powerful than, the energy loss bound. Unfortu- 
nately for rriY > 40 keV the hidden photon flux cannot be calculated from an unperturbed 
solar model and the bounds become indeterminate. 

Finally, it is interesting to compare these results with the earlier works ES]- The 
energy loss bounds presented here are weaker except in the resonant region, whose effects were 
apparently not considered before. The extremely good results of the CAST collaboration, 
are however able to reverse the situation for low masses. Even with the lower flux calculated 
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in this paper the CAST exclusion hne goes deeper, not only than the former helioscope 
bound, but also than the energy loss limit. Finally it is interesting to note that the lifetime 
bound derived here extends one order of magnitude further in mass. Therefore, except for 
two islands around m^i ~ 1 eV and m^i ~ 1 keV, the bounds presented here improve the 
earlier limits. 

Our last check should be to ensure that the 5* absorption inside the Sun is negligible as 
it has been assumed in all the above. If this were not the case, not only the hidden photon 
flux will decrease with respect to the above estimates but the same solar internal structure 
will require an unacceptable readjustment because of the resulting non-local energy transfer 
|63j . I have checked that absorption is not significant except for the case of massive hidden 
photons m^i ^ 295 eV with x ^ 10~^. The main reasons for that are that absorption in 
the low mass range is suppressed by my and the resonant regions in the intermediate mass 
regime are never too wide. For illustration purposes the line labeled "m.f.p.= Rq^ shows the 
value of the mixing parameter x for which a T-hidden photon produced in the solar center 
would have a mean free path of the order of the solar radius. This mean free path is an 
average over the radial trajectory using at every place of the Sun the highest possible value 
of the absorption rate (energies near the plasma frequency) , so it is again a very conservative 
estimate. 

In this context I should emphasize that the region above the "Lifetime" and "Helioscope" 
curves is not strictly ruled out by the energy loss and CAST limits since it assumes that 
hidden photon emission is a small perturbation of the standard solar model [37]. Only adding 
the fact that the standard solar model agrees very well with helioseismological data and the 
observed neutrino fluxes, and these are typically very sensitive to the internal structure of 
the Sun, one concludes by "reductio ad absurdum" that this region is severely excluded. 

As a final remark we should keep in mind that these bounds can be completely different 
if, in addition to the hidden photon, other low mass particles exist in the hidden sector. In 
particular, the stellar emission of particles of mass < 10 keV, charged under U{1\ does not 
vanish in the limit m^i — > [64j (actually this can be even true if these particles are much 
heavier |65]). In this case low hidden photon masses are constrained as much as masses of 
the order of the stellar temperature. 



5 Conclusions 

I have addressed the calculation of the solar emission of a hypothetical hidden sector pho- 
ton mixing kinetically with the standard model ordinary photon. I have shown that a 
resonant effect is possible when the dispersion relation of solar plasmons fits the particle-like 
dispersion relation of the hidden photon. This happens for transverse plasmons if the hidden 
photon mass m^i lies in the range 1 ~ 295 eV (the range of the plasma frequency in the 
solar model used) and for longitudinal plasmons as long as m^i < 295 eV. 
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The conservative requirement that the hidden photon luminosity should not exceed the 
solar standard luminosity bounds the amount of kinetic mixing up to x ^ 10~^^ depending 
on the mass and the polarization. At masses beyond 1 eV, where the strongest bound is 
reached, the emission of transversally polarized hidden photons dominates over the emission 
of longitudinal ones. Below this mass both polarizations contribute in a similar amount. 

At low masses the bounds are weaker, relaxing proportionally to my. However, the non 
observation of a signal in the CAST axion helioscope improves the bounds in this region up 
to 2 orders of magnitude. Altogether, these are the best limits on the mixing parameter x 
in the range 3 meV< m^i < 40 keV. A small room for improvement is available for large 
masses m^i ^ 10 keV if the CAST detectors were to rise their top energy threshold. 

It should be interesting to extend this study to other stellar objects like supernovae, white 
dwarfs, red giants and horizontal branch stars, since these can provide stronger bounds than 
the Sun, specially at masses where a resonant production is possible. 
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Appendix 



A Thomson dispersion of longitudinal plasmons 

The dispersion relation for longitudinal plasmons allows for space-like excitations. These 
plasmons can be coherently absorbed by the electrons in the plasma giving rise to an order 
0{a) absorption - the so-called Landau damping [661 [67]- which can be important for the bulk 
emission. However, only time-like L plasmons, whose "Landau" absorption is kinematically 
forbidden, can be resonantly produced making natural to restrict this discussion to them. 
The next order contribution to the absorption rate of plasmons is due to Thomson dispersion, 
7ae~ — > 7a'e~. Given that typical energies in the solar interior are much smaller than the 
electron mass, it is justified to focus on the non-relativistic limit uie —>■ oo. It is worth noting 
that in this case the transitions to transverse plasmons 7Le~ — > 7x6" are suppressed by phase 
space since uj^ik) grows much faster with k than uJi{k) and the target electrons cannot 
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transfer energy efficiently at such small energies. Therefore I only address the calculation of 
the 7Le~ — > 7Le~ process. 

Let us begin by defining the polarization vectors for the initial and final photon states. 
Writing k'^ = {uj]0,0,k) f and k'^^ = {uj'; k' sm6,0, k' cos 9) for the initial and final photon 
4-momenta we can choose 



-Xk;0,0,u) ; e 



1 



The matrix element can be written as 



- k'^ 



{k'] J sin 6*, 0, J cos( 



(39) 



M 



where 



M 



flU 



-e^^fZ^fZ' u{p', s' 



2p-k + kl ^ 2p-k' - 



u{p, s) 



(40) 



(41) 



and Z, Z' are the wave function renormalization factors of the initial and final plasmons 
given by 



Zi.{k) 



'00^ -k^ 3cu2 - 2(cu2 - t;2A;2) - P 



(41) 



where = bT/rUe <^ 1 and u, k are understood to satisfy the dispersion relation eq. (10)(L). 
In practice, because we restrict ourselves to time-like plasmons u > k, we can set Zl ~ 1. 

Because of charge conservation, M'^'^ satisfies necessarily the conditions 



k^M 



M^"^ - kM^^ = 



k'M"" 



uj'M 



k"'M' 







that can be used to express 7V1 in a very convenient way 

v/(cu2- A;2)(cu'2- A;'2) 



Mi 



[M^^ sine + A^33 



(42) 



(43) 



Note that the renormalization factors in eq. (41) will cancel the prefactor in the above 
equation. 



In the non relativistic limit we can use p^ 



P 



(m; 0) for the Dirac spinors and we 



can neglect the squared 4-momenta of the photons in the denominators of eq. (41). Now 
it is easy to evaluate A^*^. Recall that u{s')Yu{s) and are proportional to the velocity 
of the electrons and therefore are suppressed with respect to ?l(s')7°u(s) = u''{s')u{s) = 
2m^'^^. This argument is sufficient to neglect the terms proportional to p'^ in the numerators. 
Moreover, using the commutation relations {7'^,7'^} = 2g'^'^ and also (7°)^ = —(7*)^ = I5 
it turns out that terms containing three gamma matrices can be either reduced to terms 
containing only one or they appear in pairs that cancel out. Finally, Ai^^ is found to be 
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negligible compared with A4^^ c:^ 2e^^'^^ and the properly averaged squared matrix element 



IS 



= =4e^cos2^ . (43) 

s,s' 

The averaged dispersion rate is obtained by integrating the Lorentz invariant phase space of 
the final particles and averaging over the thermal distribution of electrons in the initial state, 
weighting with the appropriated stimulation/blocking factors. In practice, however, given 
that electrons are non-relativistic we can avoid their thermal average and blocking factors 
and simply multiply by the electron density rif, to get 

Tl = / . "^xf . "^xY i^^yS^P' + k'-p- k)\M\\l + n^Eiuj')) (44) 

^ 2uj2meJ (27r)32c<;'(27r)32E^^ J Ki^ ^ ;i I V ^ bev v ; 



a'^Tie f 27rdk'k'^dcos6 

LOnie J LO'lTLe 

a^Ue 27idk'k'^ ^ , Svra^rie k 

'kk 

E' 



S{E' + Lj' -m-u) cos2^(l + UBEiuj')) (45) 



f 271 dk k cos 6 , , Svra rip k , , , , 



where in the second step I have used the cos ^-dependence of the energy of the outgoing 
electron E' = + (k — k')^ = a/ + A;^ _|_ i^i"^ _ ikk' cos Q to remove the Dirac's delta 

enforcing energy conservation. Together with the dispersion relation, this establishes the 
relation 

k' 3T + cup cos Q 

T = , ' 46 

k 61 + uj-p 

that bounds the integral over k' to the small interval (fc^ = k{l — 2ti;p/3T), k). Finally note 



that the stimulation factor 1 + n-QEi!^') cancels the term in the parenthesis of eq. (11) since 
~ so we obtain 

T r 1 Ma^n^ k 
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